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I.  Abstract 

For  high  Mach  number  flows,  M  >  4,  piston  theory  has  been  used  to  calculate  the  pressures  on  the 
surfaces  of  a  vehicle.  In  a  two-dimensional  inviscid  flow,  a  perpendicular  column  of  fluid  stays  intact  as  it 
passes  over  a  solid  surface.  Thus,  the  pressure  at  the  surface  can  be  calculated  assuming  the  surface  were  a 
piston  moving  into  a  column  of  fluid.  In  this  work,  first-order  piston  theory  is  used  to  calculate  the  forces, 
moments,  and  stability  derivatives  for  longitudinal  motion  of  a  hypersonic  vehicle.  Piston  theory  predicts  a 
relationship  between  the  local  pressure  on  a  surface  and  the  normal  component  of  fluid  velocity  produced 
by  the  surface’s  motion.  The  advantage  of  piston  theory  over  other  techniques,  such  as  Prandtl-Meyer  flow, 
oblique  shock,  or  Newtonian  impact  theory,  is  that  unsteady  aerodynamic  effects  can  be  included  in  the 
model.  Prandtl-Meyer  flow  and  oblique  shock  theory  are  utilized  to  provide  flow  properties  over  the  surfaces 
of  the  vehicle.  These  flow  properties  are  used  to  determine  the  steady  forces  and  moments  and  are  also 
included  in  the  unsteady  flow  calculations.  Thus,  this  work  utilizes  a  combination  of  Prandtl-Meyer  flow, 
oblique  shock,  and  piston  theory  to  calculate  forces  and  moments.  The  unsteady  effects,  considered  in  this 
work,  include  perturbations  in  the  linear  velocities  and  angular  rates,  due  to  rigid  body  motion.  A  flexible 
vehicle  model  is  developed  to  take  into  account  the  aeroelastic  behavior  of  the  vehicle.  The  vehicle  forebody 
and  aftbody  are  modelled  as  cantilever  beams  fixed  at  the  center-of-gravity.  Piston  theory  is  used  to  account 
for  the  changes  in  the  forces  and  moments  due  to  the  flexing  of  the  vehicle.  Piston  theory  yields  an  analytical 
model  for  the  longitudinal  motion  of  the  vehicle,  thus  allowing  design  trade  studies  to  be  performed  while 
still  providing  insight  into  the  physics  of  the  problem. 

II.  Introduction 

In  the  1980’s,  the  National  Aerospace  Plane  (NASP)  program  commenced,  with  its  goal  being  a  feasibility 
study  for  a  single-stage  to  orbit  (SSTO)  vehicle,  which  was  reusable  and  could  take  off  and  land  horizontally. 
The  NASP  was  to  be  powered  by  a  supersonic  combustion  ramjet  (scramjet)  engine.  Although  this  program 
was  cancelled  in  the  1990’s,  a  great  deal  of  knowledge  was  gained  and  it  spawned  future  programs,  including 
the  hypersonic  systems  technology  program  (HySTP),  initiated  in  late  1994,  and  the  NASA  X-43A.  The 
HySTP’s  goal  was  to  transfer  the  accomplishments  of  the  NASP  program  to  a  technology  demonstration 
program.  This  program  was  cancelled  in  early  1995.  The  NASA  X-43A  set  new  world  speed  records  in  2004, 
reaching  Mach  6.8  and  Mach  9.6  on  two  separate  occasions  with  a  scramjet  engine.  These  flights  were  the 
culmination  of  NASA’s  Hyper-X  program,  with  the  objective  being  to  explore  alternatives  to  rocket  power 
for  space  access  vehicles. 
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With  renewed  interest  in  space  operations  worldwide,  there  is  an  interest  in  hypersonic  aerodynamics 
research.  The  scramjet  engine  will  likely  play  a  major  role  in  future  hypersonic  vehicles.  Unlike  a  conventional 
turbojet  engine,  a  scramjet  engine  does  not  use  high  speed  turbomachinery  to  compress  the  air  before  it 
reaches  the  combustor.  Instead,  it  relies  upon  the  rise  in  pressure  across  oblique  shock  waves  located  in  front 
of  the  inlet.  Furthermore,  the  flow  through  the  entire  engine  is  supersonic  in  contrast  to  a  ramjet,  where  the 
flow  speeds  are  subsonic  through  the  combustor.  On  configurations  like  the  NASP  and  X-43A,  the  underside 
of  the  airframe  must  function  as  the  air  inlet  mechanism  and  the  exhaust  nozzle.  Therefore,  integration  of 
the  airframe  and  engine  are  critical  to  success  of  a  scramjet  powered  vehicle. 

Scramjets  could  be  used  as  part  of  a  multi-stage  launch  vehicle  that  would  include  multiple  propulsion 
systems  to  perform  a  mission.  The  factor  driving  research  towards  scramjets  and  away  from  rockets  is  cost; 
scramjets  would  substantially  lower  costs  because  it  is  an  airbreathing  engine.  Airbreathing  engines  don’t 
require  oxidizer  to  be  carried  by  the  vehicle,  hence  increasing  the  payload  and  reducing  the  quantity  of  fuel 
carried. 

Unsteady  aerodynamics  are  a  phenomenon  that  must  be  considered  in  the  development  and  optimization 
of  future  hypersonic  vehicles.  The  combined  effects  of  a  slender  flexible  vehicle  travelling  at  high  speeds  and 
subjected  to  large  forces  may  lead  to  significant  unsteady  aerodynamic  effects.  Hence,  understanding  the 
concepts  and  consequences  of  time-dependent  aerodynamic  flows  is  critical  to  the  successful  development  of 
this  type  of  vehicle. 

Piston  theory  is  a  technique  that  has  been  used  for  years  to  determine  the  pressure  distributions  on  an 
airfoil/vehicle,  when  the  Mach  number  is  sufficiently  high.  Lighthill1  discussed  the  application  of  piston 
theory  on  oscillating  airfoils  some  50  years  ago.  Ashley  and  Zartarian2  discuss  piston  theory  while  providing 
a  number  of  examples  of  the  application  of  piston  theory  to  specific  problems.  More  recently,  Tarpley3 
discussed  the  computation  of  stability  derivatives  for  a  caret-wing  waverider  using  piston  theory,  which 
requires  the  analysis  of  unsteady  flow  over  the  vehicle.  Piston  theory  allows  the  inclusion  of  unsteady 
aerodynamic  effects  in  the  model  and  a  closed  form  solution  can  be  found  for  these  unsteady  effects. 

In  this  work,  piston  theory  is  applied  to  a  hypothetical  2-dimensional  hypersonic  vehicle  powered  with 
a  scramjet.  This  work  uses  first-order  piston  theory  to  compute  unsteady  effects  behind  shock  waves  and 
expansion  fans.  A  recent  study4  revealed  that  this  method  of  computing  unsteady  aerodynamic  effects 
delivered  highly  accurate  results  wheu  compared  to  computational  fluid  dynamics  solutions  and  higher  order 
piston  theory  models.  This  paper  builds  on  previous  work5  by  incorporating  a  flexible  vehicle  model  and 
determining  the  perturbations  to  forces  and  moments  due  to  flexibility. 

In  Section  III,  the  vehicle  analyzed  in  this  work  is  described,  while  the  steady  forces  and  pressures  on  the 
vehicle’s  surfaces  are  provided  in  Section  IV.  The  unsteady  effects  due  to  the  aerodynamic  control  surfaces 
are  included  in  Section  V,  the  total  rigid  body  forces  and  moments  are  discussed  in  Section  VI,  while  the 
aeroelastic  model  is  developed  in  Section  VII.  Results  are  provided  in  Section  VIII,  conclusions  are  given  in 
Section  IX,  and  Appendix  A  contains  detailed  calculations  for  many  of  the  flexible  stability  derivatives. 

III.  HSV  Model 

Figure  1  shows  the  2-dimensional  hypersonic  vehicle  considered  in  this  work.6  The  longitudinal  force  and 
moment  analysis  is  taken  as  unit  depth  into  the  page.  The  vehicle  consists  of  4  surfaces:  an  upper  surface 
(defined  by  points  cf)  and  three  lower  surfaces  (defined  by  points  cd,  gh,  and  ef).  All  pertinent  lengths  and 
dimensions  are  in  units  of  feet  and  degrees,  respectively.  The  total  length  of  the  vehicle  is  L  =  100 ft  and 
the  notation  for  lengths  is  Lf  =  length  of  the  forebody,  Ln  =  length  of  the  engine  nacelle,  La  =  length  of 
the  aftbody,  Le  is  the  length  of  the  elevator,  Lc  is  the  length  of  the  canard,  Xf  is  the  distance  from  the  C.G. 
to  the  front  of  the  vehicle,  xa  is  the  distance  from  the  C.G.  to  the  rear  of  the  vehicle,  xeiev  and  zeiev  are  the 
distances  from  the  C.G.  to  the  midpoint  of  the  elevator  in  the  x  and  z  directions,  respectively,  xcana rd  is  the 
distance  from  the  C.G.  to  the  midpoint  of  the  canard,  and  hi  is  the  engine  height.  The  vehicle  dimensions 
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Figure  1.  Hypersonic  Vehicle. 


are 


The  vehicle  angles  are 


L  =  100  ft 
Lf  =  47  ft 
La  =  33  ft 
Ln  =  20  ft 
Le  =  17  ft 
Lc  =  10  ft 
Xf  =  55  ft 
xa  =  45  ft 
Xelev  —  30  ft 
Zelev  =  3.5  ft 


hi  =  3.25  ft 


Xcanard  —  40  ft 


n,u  =  3° 
ti,l  =  6° 
r2  =  14.41° 


Additionally,  the  vehicle  mass  and  moment  of  inertia  are 


Mass  =  300 
Jvy  =  500, 000  slu9~ft2 


(2) 


(3) 
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and  the  mean  aerodynamic  chord  (c)  and  planform  area  ( S )  are  defined  as 


c  =  L 
S  =  L 2 


(4) 


The  goal  is  to  apply  piston  theory  to  this  2-dimensional  vehicle  geometry  to  determine  the  pressure 
distribution  on  the  surfaces  of  the  vehicle,  which,  in  turn,  can  be  used  to  evaluate  the  forces  and  moments. 
The  pressure  on  the  face  of  a  piston  moving  into  a  column  of  perfect  gas  is2 


P 

p 

1  OO 


7  ~  1  Vn  \ 

2  Uqo  / 


z-t 

7-1 


(5) 


where  the  subscript  ”oo”  refers  to  the  steady  flow  conditions  past  the  surface,  Vn  is  the  velocity  of  the  surface 
normal  to  the  steady  flow,  Ooo  is  the  freestream  speed  of  sound,  and  P  is  the  surface  pressure.  Taking  the 
binomial  expansion  of  Eq.  5,  to  first  order,  produces 

P  2y  7  -  1  Vn  7 Vn 

=  1  +  - ™=i+lPL  (6) 

MX)  7  1  2  CLqo  &OO 

Multiplying  through  by  P^  and  using  the  perfect  gas  law  (P  =  pRT)  and  the  definition  of  the  speed  of 
sound  (a2  =  jRT)  yields  the  basic  result  from  first-order  linear  piston  theory 


P  —  Poo  +  Poo  ®oo  Fi 


(7) 


where  7  is  the  ratio  of  specific  heats  and  R  is  the  gas  constant.  The  infinitesimal  force,  on  an  elemental  area, 
due  to  the  pressure  is  given  by 

dF  =  -PdAn  (8) 

where  dA  is  a  surface  element  and  n  is  the  outward  pointing  normal.  Substituting  Eq.  7  into  Eq.  8  yields 

dF  =  (-Poo  -  poodooVn)  dAn  (9) 

The  normal  velocity  can  be  computed  by  taking  the  dot  product  of  the  flow  velocity  over  a  surface  and  the 
outward  pointing  normal  for  that  surface.  Hence,  Eq.  9  becomes 

dF  =  (—Poo  —  PooO-oo  [V  •  n])  dAn  (10) 


Equation  10  is  the  basic  result  upon  which  this  work  is  based.  From  this  equation,  it  is  seen  that  in  order 
to  compute  the  forces  acting  on  a  surface,  one  must  determine  the  properties  of  the  flow  past  the  surface 
(properties  behind  a  shock,  expansion  fan,  or  freestream),  the  velocity  of  the  surface  relative  to  the  airstream, 
V,  the  outward  pointing  surface  normal,  n,  and  the  surface  element,  dA.  The  work  that  follows  will  develop 
these  quantities  for  the  upper  and  lower  surfaces  of  the  vehicle,  as  well  as  for  the  control  effectors. 


IV.  Vehicle  Surface  Pressure  Distributions  and  Forces 

The  differential  forces  on  the  surfaces  of  the  vehicle  were  computed  by  Oppenheimer  and  Dornan5  and 
are  repeated  below: 


dFcf  =  (— Pcf  —  pcfCLcf  {[n  +  5  tan  77,17  (a;  —  Xf)\  sinri^  —  [w  —  qx]  cosTpjy})  ^sinTi,cd  —  cosrpc/fcj  secTpyda; 

(11) 

dFcd  =  (-Pad  -  Podded  {[«  -  g tan rpL  ( x  -  xj)]  sinrijL  +  [w  -  qx]  costal})  sin 77, l*  +  cosn^fc  sec Ti^dx 

(12) 


dFgh  =  (- Pgh  -  Pghdgh  {w  -  qx}) 


Ik 


dx 


(13) 


dFe/  =  (- Pef  -  Pefaef  {-  ( u  +  qrefz)  sin  (ti>C/  +  r2)  +  (w  -  qx)  cos  (n, u  +  r2)})  nef  sec  (77, j/  +  r2)  dx 


(14) 


4  of  25 


American  Institute  of  Aeronautics  and  Astronautics 


where 


(15) 


ef  =  -  sin  (ri,c/  +  T2)i  +  cos  (t^u  +  r2)  k 


and 


refz  =  [tan  (t1}U  +  r2)  {x  +  xa)  -  Ltanr1>c/] 


(16) 


where  u,  w,  and  q  are  small  perturbations  in  a  steady  flight  condition  of  the  forward  velocity,  normal  velocity, 
and  pitch  rate,  respectively.  The  steady  forces  are  computed  by  integrating  the  steady  components  of  the 
differential  forces  over  the  corresponding  surface.  Performing  the  integrations  yields 


F cfa  =  -PcfXa  sec  ThU  sinryj/z  -  cos  ntuk 
F cff  =  —  PcfXf  secTitu  sinri^z  —  cosri^fc 


Xcfai  fa  k 

Xcffi  T  Zcff  k 


(17) 


where  Xcfa ,  Zcfa  are  the  components  of  the  aftbody  upper  surface  force  in  the  x  and  z  directions,  respectively, 
and  Xcff ,  Zcff  are  the  components  of  the  forebody  upper  surface  force  in  the  x  and  z  directions.  The  lower 
surface  forces  are 


Fcd  =  -PcdLf  sec 


sin  T\  L,i  +  cos  t\  Lk 


—  Xcdi  P  Zcdk 


F  gh  —  PghLnk  —  Xgfoi  -\-  Zgfok 


(18) 


while  on  the  rear  ramp,  the  vector  force  due  to  the  external  nozzle  is 


F  ef  = 


Ln.PePry. 


In 


Pe 


COs(n,C/  +  T2)  (Pe  -  Poo)  Pc 


sin(nic/ +  r2)  z  -  cos  (tiiC/ +  r2)  fc  =Xefi  +  Zefk  (19) 


where  Xef  and  Zef  are  the  axial  and  normal  force  components  of  the  external  nozzle  force  and  Pe  is  the 
engine  exit  pressure. 


A.  Control  Surfaces 

The  control  surfaces  are  an  elevator  and  canard  as  shown  in  Fig.  1.  Both  surfaces  are  modelled  as  flat  plates 
hinged  at  their  midpoints  so  the  entire  surface  deflects.  Positive  deflection  is  defined  as  trailing  edge  down. 
The  x  and  z  positions  of  the  midpoint  of  the  elevator  and  canard  referenced  to  the  C.G.  are  xeiev,  zeiev  and 
x canard:  0.  The  canard  was  added  to  this  vehicle  for  two  reasons.  First,  the  elevator  produces  a  significant 
amount  of  lift  and  results  in  a  non-minimum  phase  flight  path  angle  response.'  Adding  the  canard  and 
forcing  it  to  respond  in  conjunction  with  the  elevator  compensates  for  the  lift  produced  by  the  elevator  and 
reduces  the  non-minimum  phase  behavior.  Secondly,  the  canard  can  be  utilized  for  low  speed  control  and 
reduces  the  takeoff  and  landing  speeds  since  they  produce  positive  lift  in  conjunction  with  a  nose  up  pitching 
moment.  In  this  work,  the  canard  is  not  ganged  to  the  elevator,  rather,  it  is  free  to  move.  Hence,  the 
non-minimum  phase  flight  path  angle  response  is  still  present.  However,  this  allows  the  determination  of  the 
steady  and  unsteady  forces  and  moments  due  to  the  canard. 


1.  Elevator 

The  upper  and  lower  forces  on  the  elevator  can  be  computed  as 


1*  3'elev~^~  <2  COS  <5e 

F eu  =  1  [~ Peu  ~  Peuaeu  {^ecr  ■  neiy}] 

—  sindez  —  cos  5ek 

J  X elev  COS  <5e 

sec  5edx 


r—Xelev+kf  COS  <5e 


F  — 

r  eL  — 


'  ~ Xeiev - #  COS  5e 


[  Fei  peLaeL  {XeL  .nei}] 


sin<5ez  +  cos  6ek 


sec  Sedx 


where 


Ve[/  .  neu  =  -(u-  q  {zeiev  +  tan  5e  {x  +  xeiev)})  sinde  -  (w  -  qx)  cos  6e 
Vei  .  n eL=(u-  q  {zeiev  +  tan  Se  (x  +  xeiel^)})  sin  5e  +  (w  -  qx)  cos  Se 


(20) 

(21) 


(22) 

(23) 
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2.  Canard 


In  order  to  determine  the  steady  and  unsteady  forces  and  moments  due  to  the  canard,  Eq.  10  must  be 
evaluated.  Hence,  the  flow  properties  on  the  top  and  bottom  surfaces  of  the  canard  along  with  a  position 
vector,  normal  vector,  and  surface  element  must  be  determined.  For  the  canard,  the  flow  properties  are 


If  Sc  =  — a  =>■  Freestream 
If  sc  <  —a  =>  Shock  on  top,  expansion  on  bottom 
If  Sc  >  —a  =>■  Expansion  on  top,  shock  on  bottom 

The  position  vector  from  the  C.G.  to  an  arbitrary  point  on  the  canard  is  given  by 

Yc  —  xi  A  [  tan  Sc  {x  Xcanar(i )  +  Z canard ] 

%  canard  f  COS  Sc  ft  %  C  % canard  T  f  COS  Sc 


(24) 


(25) 


where  zcanard  =  0  for  this  configuration.  The  outward  pointing  normal  vectors  for  this  surface,  relative  to 
the  body  axis,  are 


and  the  surface  element  is 


ncu  =  —  sin  Sci  —  cos  Sck 
nCr  =  sin  6ci  +  cos  5ck 


dAc  =  sec(5c(l)da; 


(26) 


(27) 


To  compute  the  forces,  moments,  and  stability  derivatives,  consider  small  perturbations,  from  a  steady  flight 
condition,  in  the  velocities  u  and  w  and  the  pitch  rate  q.  The  velocity  of  a  point  on  the  upper  and  lower 
surfaces  of  the  canard,  due  to  these  perturbations,  is 


VC[/  =  ( Vcu  cos  Sc  +  u)  i  +  (— Vcu  sin  Sc  +  w)  k  +  u>  x  rc 
VCL  =  ( VCL  cos  6C  +u)i  +  (— VCL  sin  Sc  +  w)  k  +  uj  x  rc 


(28) 


where  Vcu  and  VCL  are  the  flow  velocities  on  the  upper  and  lower  surface  of  the  canard  (as  determined  from 
the  flow  analysis)  and  uj  =  qj  is  the  angular  rate  vector.  Using  Eqs.  24  -  28  in  Eq.  10  yields  the  differential 
forces  due  to  the  canard: 


dF cu  —  [  Pcu  Pcu^cu  {  sin 5C  ( u  T  q  [  tandc  ( x  xcanar(f)  -\-  zcanar(f\)  cosSc  (w  (72)}]  xvcudAcu 
dFCL  —  [  Pcl  PcLaCL  {sin 5C  (u  T  q  [  tan<5c  ( x  xcanar(i)  T  Zcanardl )  T  C0S(5c  ( w  mCLdACL 


The  steady  forces  are  computed  by  integrating  the  steady  component  of  Eq.  29,  hence, 


(29) 


F-  =  l 
fcl  =  r 


'^canard~^~  2°  ^ c 

^canard - ^  COS  Sc 

X canard~\  ^  COS  5C 


—Pen  (  _  sin  5ci  —  cos  Sck  )  sec  6cdx  =  PCIT Lc  sin  Sci  +  PCTt Lc  cos  5ck  =  XCTTi  +  ZCTT  k 


Lc  _  x  -PcL  ( sin  5ci  +  cos  Sck )  sec  5cdx  =  —PCL  Lc  sin  Sci  —  PCL  Lc  cos  Sck  =  XCL  i  +  ZCL  k 


Xcanard  o°  COS  8, 


The  corresponding  moments  are  calculated  using  MCu  =  rc  x  FCu  and  MCl  =  rc  x  FCl  such  that 

Mcu  %canardP 2[/TcCOS$c  T  Z canar d P c\j  Zc  sin  Sc  —  X canard P c\j  COS  6 c 
MCl  =  ^  canard  Pcl  c  COS  Sc  Zcanar^PCLPc  sin  Sc  =  Xcanar^PCLPc  COS  5C 

where  the  last  inequalities  in  Eq.  31  result  from  zcanard  being  zero. 


(30) 


(31) 


V.  Unsteady  Control  Surface  Effects 

The  unsteady  effects  due  to  velocity  and  rate  perturbations  on  the  control  surfaces  are  found  using  the 
unsteady  portions  of  Eqs.  20,  21,  and  29. 
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A.  Elevator  Unsteady  Effects 

The  normal  force  stability  derivative  due  to  the  elevator  is 

(CzJs -  =  ^s  I  {dF)z~w  =  ^S  I  {dFeu)*~™  +  ^sf  (dFei)*-« 


Performing  the  necessary  substitutions  produces 


(CzJSe  =  — g 


f*  elev~\~  "2  COS  8e 


Qoo  S  I  J -xeiev  -  Af  cos  8e 


—peu  Cleu  W  COS  5edx  + 


-xeiev  +  :^-  cos  <5e 


‘Eelev  2  C°S 


— peLaeLw  cos  Sedx  (33) 


Integrating,  letting  w  ~  V^a,  and  simplifying  yields 


dcz\  -(p 

euaeu  d”  PeLaeL  )  Uoo  Le  COS2  Se 


Using  similar  analysis,  the  axial  force  component  becomes 


{CXw)s*  “  goo S  /  {dF)*~w  ~  goo5  /  +  9oo5  /  (dFeJ— 

&C x  \  _  (<Peu  d"  Pe l  Q'e l  )  ^oo  -^e  COS  Se  sin  Se 


The  pitching  moment  contribution  due  to  w  motion  is 


(CMJSe  =  [/  z(dFeu)x_w  +  J  z(dFeL)x_w-  J  x(dFeu)z_w-  j  x(dFeL)z_w 

{ dC]\/[\  +  peLaeL )  {^zeievLe  sin  Se  cos  Se  xeievLe  cos  $e)  Voo 


The  normal  force  increment  due  to  pitch  rate  is 

(Cz,),.  =  ^5  [/ 

,  ^  DC z  \  (Pe(/^e[/  T  [Le  sin  Se  COS  (5e  T  Xelev  tan  (ie)  a^eiei.ie] 

V  ^g  Ae  fe5. 

The  final  stability  derivative  for  the  elevator  is  the  pitching  moment  increment  due  to  pitch  rate,  which 
becomes 

(Cm, )se  =  ~  ^  J  z  (dFeu)x_q  +  J  z  (dFeL)xq  —  J  x  {dFeu)z_q  —  J  x  (dFeL)z_q  (38) 

Performing  the  required  operations  yields 


dCM\ 

)s. 


( Peuaejj  +  Pe_Laei) 


[Le  sin  5e  cos  Se  (tan Se  {x2lev  -  z2lev}  +  xeievzeiev  {l  -  tan2  <5e })] 


( Peu  ^eu  "b  Pe_£,®ei) 


+  d — - — u  s_  L — —  ( zeiev  +  xeiev  tan  Se)  xeievLe  tan  8e  (39) 


{Peu^eu  Pcl^gl)  ( tan  Se  1 
qoo  Sc  \  3  cos  A 


%elev  H-  2  COS  Se  /  \  %elev  2  COS  (5e 


B.  Canard  Unsteady  Effects 

The  normal  force  stability  derivative  due  to  the  canard  is 


(Cfc)*- -ss /<"«')”• +  ss /<"“)• 
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which  becomes 

{  dC  z  \  (PcjjQ'cu  PclO'cl)  (V^LcCOS 

\  da  )  Sc  q0 0S 

The  axial  force  increment  due  to  motion  in  the  vertical  direction  becomes 

(Cxw)Sc  =  J  idF cv)x-w  +  J  ( dFCL)x_w 

^ *  {  dC x  \  (pcv Q>cu  +  Pcl  Q’cl  )  VoqLq  cos  Sc  sin  5C 

V  da  )Sc  qooS 

The  pitching  moment  contribution  due  to  w  motion  is 

(^M„)s  = 


q^Sc 


J  z  (dFcu )x_w  +  J  z  [dF CL) x_w  Jx(dFcu)z_w  J  x  (dFCL)z_w 

f  3Cm\  ( PcijQ'Cu  +  Pcl^cl)  (  ^canardJ^c  Sill  Sc  COS  Sc  +  XcanardLc  COS  (5c)  Voc 

V  da  )5  <7oo  Sc 


The  normal  force  increment  due  to  pitch  rate  is 


(&,)<  = 


qooS 


j  ( dFcu)z_q  +  j  (dFCL)2 


(  dC*z\  {Pcu^cjj  +  Pcl^’Cl')  \XcanardFc  COS“  Sc  ZcanardLc  sin  Sc  COS  Sc 

V  dq  J  s  q^oS 


(41) 


(42) 


(43) 


(44) 


The  final  stability  derivative  for  the  canard  is  the  pitching  moment  increment  due  to  pitch  rate,  which 
becomes 


(cwj,  = 


i 


5c  qooSc 

Performing  the  required  operations  yields 


J  z  (dFCu)x_q  +  J  z  (dFCL)x_q  jx{dFCu)z_q  J  x(dFCL)z 


(45) 


(' °CuacU  +  PcLacL  ) 

sin^  SCL ^  2  2  2  2  2  2  cos^ 

V  dq  J 

Sc  QooSc 

12  zcanardLc  sm  sc  ^  sin  5C  cos  Sc  xcanardLc  cos  Sc  ^ 

(46) 


VI.  Total  Forces  and  Moments  -  Rigid  Body 

With  the  inclusion  of  the  stability  derivatives,  the  thrust,  and  resulting  engine  moment,5  the  total 
aerodynamic  forces  and  moments  on  the  vehicle  are 


Xtotai  —  Xcff  +  Xcfa  +  Xcd  +  Xgh  +  Xef  +  XeL  +  Xeu  +  Xcu  +  XCL  +  T  +  qooS 


dCx 


da 


a 


Ztotal  —  Zcf;  +  Zcfa  +  Zed,  +  Zgh  +  Zef  +  ZeL  +  Zeu  +  ZCJJ  +  ZCL  +  qooS  ^  a  +  qooS 

,  c(dCz\  ,  c(dCz\  ,  a(dCz\  qc  ,  JdCz\  qc 

+q-s  {- da)se  a + q~s  Ur>L a + q-s  y^r)  + q°°s  y^r)  Sc 


(47) 


(48) 


riC1 

Mtotai  =  Mcfs  +  Mcfa  +  Mcd  +  Mgh  +  Mef  +  MeL  +  Meu  +  Mcu  +  MCL  +  Mengine  +  qooSc  ^  a+ 


c.dCM  qc  ,  c_  (dCM\  ,  c.(dCM\ 
qooSc — — ——+900601— —  a  +  qooSc  — —  a 

dq  2Voo  V  da  Js  \  da  J  Sc 


a-(dCM\  qc  ,  c.(dCM\ 

q°°‘ Sc\~w),  ScUrA 


qc 


(49) 


where  Me 


and  T  are  the  moment  and  thrust  produced  by  the  engine.  It  can  be  seen  that  Eqs.  47-  49 


do  not  contain  any  inlet  turning  forces/moments  that  was  present  in  previous  work.  Here,  a  reflected  shock 
from  the  engine  inlet  was  modelled  and  performs  the  flow  turning. 
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VII.  Flexible  Effects 


Thus  far,  only  unsteady  effects  due  to  rigid  body  motion  have  been  considered.  Full  scale  hypersonic 
airbreathing  vehicles  are  expected  to  be  long  and  slender  and  thus  highly  flexible.  This  structural  bending 
affects  downstream  flow  resulting  in  localized  changes  in  surface  pressure  along  the  body  and  thus,  should  be 
incorporated  into  the  model.  In  this  section,  the  flexible  effects  are  included  in  the  analysis.  Piston  theory 
is  still  used  to  determine  the  pressure  distribution  on  the  surfaces  of  the  vehicle  and  much  of  the  analysis 
already  presented  can  be  easily  adapted  to  included  these  additional  effects. 


A.  Flexible  Model 


In  order  to  develop  the  aeroelastic  model,  a  few  assumptions  are  made.  First,  the  flexible  vehicle  is  modelled 
as  two  cantilever  beams  fixed  at  the  c.g.  (one  for  the  forebody  section  of  the  vehicle  and  one  for  the  aftbody 
section  of  the  vehicle).  Second,  the  beams  are  assumed  to  have  constant  mass  density,  area,  and  flexural 
rigidity  (El),  where  El  is  chosen  to  give  the  desired  natural  frequency  of  vibration.  Also,  it  is  assumed 
that  the  flexible  effects  only  perturb  the  surface  velocities  in  the  z  (normal)  direction.  This  assumption  is 
justified  using  the  small  angle  approximation,  i.e. ,  the  deflection  of  the  tip  of  the  beam  is  small  compared 
to  the  length  of  the  beam.  Lastly,  it  is  assumed  that  the  change  in  angle  of  attack,  as  seen  by  the  entire 
forebody,  is  the  change  in  angle  of  attack  as  seen  by  the  tip  of  the  vehicle,  point  (xf,  0)  in  Fig.  1.  Since  the 
tip  experiences  the  largest  deflection,  this  is  a  worst-case  assumption.  This  change  in  angle  of  attack  is  used 
to  compute  the  flow  properties  behind  the  bow  shock  and  those  flow  properties  are  assumed  constant  over 
the  lower  forebody. 

The  transverse  vibrations  in  the  beam  satisfy  the  following  partial  differential  equation:8 


T^Td4w(x,t)  „  d2w(x,t) 


=  0 


(50) 


where  w(x,  t)  describes  the  position  of  the  beam,  relative  to  the  body  x-axis,  E  is  Young’s  Modulus,  I  is  the 
moment  of  inertia  of  the  beam  cross-section  about  the  y-axis,  and  rh  is  the  mass  density  of  the  beam.  This 
problem  is  typically  solved  using  separation  of  variables.  Assume 


w{x,t)  =  $(x)r](t) 


(51) 


Substituting  the  expression  for  w(x,t)  in  Eq.  51  into  Eq.  50  and  simplifying  yields 

El  «94<&(rr)  = _ 1_  d2rj(t) 

rh${x)  dx4  r/(t)  dt 2 

Since  the  left  side  of  Eq.  52  does  not  change  as  time  varies,  the  right  side  of  Eq.  52  must  be  a  constant. 
Similarly,  since  the  right  side  of  Eq.  52  does  not  change  as  x  varies,  the  left  side  of  Eq.  52  must  be  a  constant. 
Let  this  constant  be  w2,  such  that 

El  d4<j>Qr)  _  1  d2V(t)  _  2 

m$(x)  dx4  T){t)  dt2 


Utilizing  the  method  of  separation  of  variables,  Eq.  53  can  be  written  as  two  differential  equations,  one  with 
respect  to  position  and  one  with  respect  to  time: 


d4$(x) 

dx4 


fi4$(x)  =  0 


where  (34 


uj2rh 

El  ‘ 


d2rj(t) 
dt 2 


+  U}2r/(t)  =  0 


The  general  solution  to  Eq.  54  is8 


(54) 

(55) 


$(x)  =  Ci  sin  /3x  +  C2  cos  fix  +  C3  sinh  fix  +  C4  cosh  fix 


(56) 


Using  the  following  boundary  conditions  for  the  forward  beam 


d2&(x)  I 
dx2  I 


$0)1 


|  X=X f  —  0 
=  0 


d3&(x) 
dx 3 
d<&(x) 
dx 


|  X=X f  —  0 

|x=0  =  0 


(57) 
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which  state  that  the  bending  moment  and  shear  force  are  zero  at  the  free  location  ( x  =  Xf)  and  the 
displacement  and  slope  are  zero  at  the  fixed  location  ( x  =  0),  along  with  the  modal  shape  expression, 
(Eq.  56),  and  simplifying  results  in  the  frequency  equation 

cos  (3fXf  cosh /3f5if  =  —  1  (58) 

Eq.  58  has  an  infinite  number  of  solutions  with  the  first  few  given  by 


PftrXf  =  1.8751, 4.6941, 7.8548, 10.9955, 14.1372, . . . 


(59) 


The  values  of  f3fjT,  r  =  1,  2, 3, . . .  in  Eq.  59  are  called  the  eigenvalues.  Corresponding  to  these  eigenvalues, 
the  natural  modes  of  the  forward  beam  are8 


$ f,r{x )  =  Af<r  [(sin  j3f>rXf  —  smhf3f,rXf)  (sin  j3f>rx  —  sinh/3/>ra;)] 
+AftJ.  [(cos (3f,rXf  +cosh f3f,rXf)  (cos f3fTx  —  cosh/l^cr)] 

where  Af^r  is  a  normalizing  factor,  selected  such  that 

/  r(x)dx  =  1 

Jo 

where  rh  f  is  the  mass  density  of  the  forebody  beam  defined  as 


Thus,  Af<r  becomes 


Af,r  — 


fa  =  Mass  f1  -  if) 

f  Xf 

1 


Af,rPl  +  ASxp2  +  Af,r 


where 


a  (sin /3f  rXf  —  sinhfif  rXf)2  r 

Af^  ~ - 4^ - [Mlf] 


(sin  f3 f  rXf  —  sinh  (3 f  rXf)  (cos  (3 f  rXf  +  cosh  (3 f  rXf)  .  . 

Af>rp2  =  2  2  w,  2  2  (sn 


■  (sin  (3frXf  —  sinh  ftfrXf)2 

Pf,r 

(cos  /3ftrXf  +  cosh  /3ftrXf)2 

xi  f.r  ^  m 


4/?/,  r 


[M3/] 


and 


(60) 

(61) 

(62) 

(63) 

(64) 

(65) 

(66) 


Mi f  =  —  2  cos  (3f,rXk  sin  (3 ftTXk  +  sinh  2 j3f<rXk  —  4  sin  j3f>rXk  cosh  /3f>rXk  +  4  cos  (3f,rXk  sinh  /3ftrXk  (67) 

M$f  =  2  cos  (3ftrXk  sin  (3f,rXk  +  sinh  2/3f,rXk  —  4  cos  j3f>rXk  sinh  f3ftrXk  —  4sin/3/iT.Xfc  cosh  (3f,rXk  +4,/3f,rXk  (68) 
For  the  aft  beam,  the  boundary  conditions  are 


d2&(x)  | 
dX  2  \x=x, 


$(aOU=o 


d3$(x) 
dx 3 
d&(x) 
dx 


|af=xa  ^ 
|x=0  =  0 


(69) 


and  the  frequency  equation  becomes 

cos  f3axa  cosh  f3axa  =  -1 


with  solutions 

f3a,rxa  =  1.8751,4.6941,  7.8548, 10.9955, 14.1372, . . . 
The  natural  modes  of  the  aft  beam  are 


(70) 

(71) 


$a,r(x)  =  Aa  r  [(sin f3airxa  -  sinh/30jrx0)  (sin/laira;  -  sinh/la>ra;)] 
3~Aa^r  [(cos  (3a  rxa  cosh  f3a,rXa')  (cos  (3a^rx  cosh/la  rx)] 


(72) 
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(73) 


and  Aa  r  is  selected  such  that 


and 


Thus,  Aair  becomes 


where 


ma$lr(x)dx  =  1 
Mass  (l  — 


A 


An 


Aa,rPl  T  Aa^rp^  A  Aa  r f 


(sin (3a,rXa  -  sinh/3aira/) 


[Mi 


A„ 


'  P2 


(sin  l3a,rXa  sinh/3a  r:ra)  (cos  fda  rxa  A  cosh/?a  rxa)  ,  .  _  -to  —  \2 

- (smpa>ra;a  -  smh  pa,rxa) 


A, 


a,rp 


(Oa,r 

(COS  {3a,rXa  A  COSh  0<irX  n  ( 

4  (3a, r 


[M3a 


(74) 

(75) 

(76) 

(77) 

(78) 


and 


M\a  2  cos  Pa_rxa  sin  Oqj-Xq  A  sinh  2f3a^rxa  4  sin  f3a^rxa  cosh  (da^rxa  -b  4  cos  Pa,rxa  sinh  (da,rxa  (79) 
M3a  2  COS  0a,rXa  Sin  {3a,rXa  A  sinh  2f3a,rXa  4  COS  Pa,rXa  sinh  (3a,rXa  4  sin  f3a,rXa  COSh  fia,r%a  A  ^(3a,rXa  ($0) 

B.  Forced  Response 

Let  the  forcing  function  in  Eq.  50  consist  of  distributed  and  concentrated  loads  so  that  Eq.  50  can  be  written 
as 

04w(x,t)  „  d2w(x,t)  ,,  x  „  ,  x  x 

;  =f(x,t)  +  Fj(t)S{x-xj) 


El¬ 


m- 


Ox 4  '  ""  dt2 

where  6(x)  is  the  dirac  delta  function  defined  as 


1  if  x  =  Xj 
0  if  x^  xj 


5(x~Xj)  = 

From  the  expansion  theorem,  the  solution  to  Eq.  81  is 

OO 

Wf(x,  t)=^2  ^/,r(ar)ry/,r(i) 

r=  1 

for  the  forebody  beam  and 

OO 

Wa(x,  t)  =  ^  &a,r(x)'lla,r(t) 


(81) 


(82) 


(83) 


(84) 


for  the  aftbody  beam  where  i]f,r(t),  rjair(t)  are  the  generalized  modal  coordinates,  for  the  forebody  and 
aftbody  beams,  that  satisfy 


Vf.r (f )  T  A  0) f ^V) f ,r(t)  —  Nf^r(t) 

Va,r(t)  A  2<(^arCJarfjar(t)  -\-  0Ja  rTJa  r(t')  —  Na  r(t ) 


(85) 


Here,  Nf  r(t),  Na  r(t)  are  generalized  forces  for  the  rth  mode  shape  of  the  forebody/aftbody  beam,  defined 
by8 


rx  f  n 

Nf,r(ty—  /  *&f,r(x)ff(x,t)dx  A  ^  ^*&f,riXj)Ff,j(t) 
J0  ,  =  1 
rxa  n 

Na,r(t)  =  /  ^a,r(x)fa(x,t)dx  +  y2<Pa,r(xj)Fa,j(t) 
J0  J=1 


(86) 
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where  n  is  the  number  of  concentrated  loads  on  the  beam.  Given  the  loading  on  the  forebody  and  aftbody 
beams,  the  generalized  forces  for  the  first  mode  become 


rxr  rxf 

A(/,  1  (t)  —  /  Qfi(x)Pcfdx  /  d*/,/ (x) Pcddx  T  4*/  i (xcanard)  ( PCjj  Pcl)Lc 

J 0  J x f—Lf 

pxa  rxa 

Na,i(t)=  $atl(x)Pcfdx-  $a}1(x)Pefdx  +  <i>atl(xeiev)(Peu  -  PeL)Le 

Jo  Jxa—La 


(87) 


These  expressions  are  displayed  in  Appendix  A. 

In  order  to  incorporate  aeroelastic  effects  into  the  model,  a  few  simplifying  assumptions  are  made.  First, 
the  vehicle  does  not  stretch  or  compress  along  the  x-axis.  Second,  for  small  displacements,  when  the  beams 
(vehicle)  flexes,  there  is  no  change  in  the  x  direction  displacement.  With  these  assumptions,  aeroelastic  effects 
only  occur  in  the  z-direction.  Additionally,  it  is  assumed  that  the  engine  nacelle  is  rigid.  The  aeroelastic 
effects  can  be  accounted  for  by  taking  the  time  derivative  of  Eqs.  83  and  84 


OO 

Wf{x,  t)  =  Y^  ®fAx)VfAt) 

r=  1 
oo 

Wa(x,t)  =  y^i^a,r{x)fla,r(t) 
r=  1 


(88) 


and  including  this  effect  in  the  expressions  for  the  velocities  on  the  upper  and  lower  surfaces,  namely  Eqs.  11, 
12,  and  14.  The  differential  forces  on  the  surfaces  become 


dFcff  =  ( —Pcf  ~  Pcfacf  {[w  +  gtanryc/  ( x  -  Xf )]  sinriiC/  -  [w  -  qx  +  Wf(x,t)]  cosryt/})  secriiC/nc/(ia: 

(89) 

dFcfa  =  (~Pcf  ~  Pcf cicf  {[«  +  gtanri;C/  (x  -  Xf)]  sinryt/  -  [w  -  qx  +  wa(x,t)\  costi>C/})  sec  TitUncfdx 

(90) 

dFcd  =  (- Pcd  -  PcdCicd  { [it  -  9  tarn'll  (x  -  xf)]  sin  t1>l  +  [w  -  qx  +  wf(x,  <)]  cosrliL})  secTf^n^dx  (91) 

dFe/  =  (~Pef  ~  Pefdef  {-  (u  +  qrefz)  sin  (rqy  +  r2)  +  (w  -  qx  +  wa(x,  t))  cos  (ri;C/  +  r2)})  sec  (t^v  +  r2)  nefdx 

(92) 

where  the  upper  surface  force  has  been  split  into  forebody  and  aftbody  parts  to  account  for  the  two  beam 
structural  model.  With  these  differential  forces,  stability  derivatives  due  to  the  bending  of  the  vehicle  can 
be  determined.  For  the  normal  force  on  the  forebody  beam, 


(Cz) 

1 

Qoo  S  , 


W  f  orebody 


qoo  S 


rXf 

fXf 

/  (^0/)^  + 

/  (dF  cd)z-w 

Jo 

>Xf-Lf 

Xf  cXf 

— Pcf  Clef  Ul  COS  T\jjdx  +  I  -pcdClcdW  cos  Ti^dx 


(93) 


C/-L/ 


Substituting  Eq.  88  into  Eq.  93  produces 

1 


(CzA 


W  forebody  q^S 


fXf  ^  f‘X  f  _  _ 

/  Pcf  Clef  E  $ f,r(x)r)f,r(t)  cosTitUdx  +  /  -podded  ^  ®  f  Ax)v f  AO  cosrltLdx 

J0  r- 1  Jxf—Lf  r=1 

(94) 


which,  for  the  first  bending  mode  becomes 

*/ 


dCz 

dry 


Cz  ...  1  r 

7/,i  <looS  Jo 


-pcfdcfQf, i(x)  cos  n, udx  + 

1 


J- 

J  X  f 


—Pcddcd^f,  i(x)  cos  n,  Ldx 


qooS 


— '2pcfacf  cosTqc/A/  i  .  .  _  .  ,  a  _  , 

— —  (sm/3/,ix/  —  smh/3/,ix/) 


A: 


Pcddcd  cos  ti,l  A/,i  ^sjn^ia,^  _  sinh  (3f,iXf)  (cos/3/,i  (x/  —  Lf)  +cosh/3/,i  (x/  —  Lf))] 


(95) 


/3f,iqooS 

PcdCicd  COSTl, lA/,i 

/3/,  i  qoo  S 


[(cos +  cosh/3/,ix/)  (—  sin/3/,i  (x/  —  Lf)  +  sinh /3/,i  (x/  —  Lf))] 
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where  A  f  i  is  the  normalizing  factor  associated  with  the  hrst  mode.  Using  the  assumption  that  the  flexible 
effects  only  perturb  the  surface  velocities  in  the  z-direction,  the  axial  force  stability  derivative  associated 
with  the  flexible  effect  becomes 

dCx 


<%,  i 


=  0 


forebody 


For  the  pitching  moment, 


(96) 


^M)w  forebody  q^Sc 


rxf  rxf 

/  z(dFcf)x-yJ+  z  (dFcd)x_ 

JO  J  Xf  —  Lf 


qooSc 


/  X  (dFcf)z-w  X  (dFcd)z_yj 

/  0  J x f  —  Lf 


(97) 


For  the  first  mode  only,  this  expression  is  displayed  in  Appendix  A. 
For  the  aftbody  beam,  the  Z  direction  force  due  to  flexibility  is 


(cz)^ 


1 


lb aftbody 


Qoo  S 


IJ  0 


(^c/),_„  + 


'  xa—La 


(^/)2_, 


and  this  expression  is  evaluated  in  Appendix  A. 

The  pitching  moment  coefficient  due  to  aftbody  flexibility  is 


(98) 


(CmU 


1 


w  aftbody  q^Sc 

1 

qooSc 


uo 


z(dFcf)x_^  + 


’(dFcf)z^  + 


,-L. 


/  xa—La 


2  (dFef)x-w 

x(dFef)z_^ 


(99) 


and  this  expression  is  evaluated  in  Appendix  A  for  the  first  bending  mode. 

With  the  inclusion  of  the  flexible  stability  derivatives,  the  total  aerodynamic  forces  and  moments  on  the 
vehicle  are 


Xtotai  —  Xcff  +  Xcfa  +  Xcd  +  Xgh  +  Xef  +  XeL  +  Xeu  +  XCu  +  XCL  +  T  +  qooS 


dCx 

da 


-a 


Fqoc  ^ 


dCx 

da 


a  +  qooS 


, dCz 


dCx 

da 


(100) 


Zt.otal  —  Zcft  +  Zcfa  +  Zcd  +  Zgh  +  Zef  +  ZeL  +  Zeu  +  Zcu  +  ZCL  +  q^S  a  +  qooS yy 


T^oc  $ 


dCz\ 

da  ). 


a  +  qooS 


dCz\  ,  JdCz 

s^),  °+fcS(^y 


qc 

5e 


Mtotai  —  Mcff  +  Mcfa  +  Mcd  +  Mgh  +  Mef  +  M£l  +  Meu  +  Mcu  +  MCL  +  M , 
_dCM  qc 


,  ..  0dCz  ,  0dCz 

+9ooA  o.  .  qood  ~ . 

dr/a,  l  oijf!  i 

,  Q-dCM 
+  qoobc— — a+ 


engine  ■  itxj  v"  Qq 

dCM  \  qc 


c-dCM  qc  ,  G.(dCM\  ,  a-(dCM\  ,  a.fdCM\  qc  ,  a.(dCM\  qc 

5“Sc_ar  ^ +  "“Sc  (■ *r) . “  ■ +  "“Sc  (■ ar) . “  +  ,“Sc  (-8T J .  wz  +  s“Sc  {sr ).  m 


c.dCM  .  ,  Q-dCM  . 

+qocdc— — r?Q>  1  +  qoodc— — ?7/,i 
UT)a,l  UWft  1 

(102) 


C.  Equations  of  Motion  for  a  Flexible  Vehicle 

The  equations  of  motion  for  the  flexible  vehicle  are6 
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a  = 


where 


Vt  =  ^  (T  cos  a  —  D)  —  g  sin  (9  —  a) 

^  (-T  sina  -  L)  +  Q  +  ^  cos  (0  -  a) 

IyyQ  -  'i  ff'lf  -  i’afla  =  M 

e  =  Q 


kfVf  +  2 CfUfVf  +  v2fVf  =  Nf  - 

kaija  +  KaU)aTla  +  ^Va  =  Na  ~ 

h  =  Vt  sin  ( 6  —  a) 


'bf^giia 

\I>  f  'arj f 
Jyy 


=  fof  x®fAx)dx 
kf  =  i  +  p- 

J  _  lyy 

=  fo  “  X$a,l(x)dx 

ka  =  1  + 


(103) 


(104) 


Clearly,  the  flexible  effects  are  coupled  into  the  pitch  rate  equation.  In  addition  to  this,  the  bending  of  the 
structure  has  an  effect  on  the  angle  of  attack  of  the  vehicle.  Since  engine  performance  is  a  function  of  shock 
angle  and  shock  angle  is  a  function  of  angle  of  attack,  a  significant  change  to  the  vehicle’s  performance  can 
occur  due  to  structural  bending.  It  is  assumed,  from  the  point  of  view  of  the  bow  shock,  that  the  entire 
forebody  observes  the  same  change  in  angle  of  attack  as  seen  at  the  nose  of  the  vehicle.  In  other  words, 
the  worst  case  change  in  angle  of  attack  is  used  for  the  entire  forebody.  This  change  in  angle  of  attack  is 
computed  as: 

Aa  =  arctan  [<!>)■  (x/)  T)f,i(t)]  (105) 


The  control  surfaces  also  see  a  change  in  angle  of  attack,  which  affects  the  properties  of  the  flow  over  the 
surfaces.  These  changes  are  given  by 


AamnBr(i  =  arctan  [<fy  ( xcanard ) 

A«cs  =  arctan  [$(>  ( xeiev )  rya,i(t)] 
so  that  the  total  angle  of  attack  seen  by  the  control  surfaces  is 

(X canard,  CT  T  A(Xcanard 

&elev  =  (X  T  A(Xeiev 


(106) 

(107) 


(108) 


The  forces  due  to  the  control  effectors,  as  given  in  Eqs.  20,  21  and  30  are  written  in  the  rigid  vehicle’s 
body  axis  frame.  However,  when  the  structure  bends,  these  forces  are  not  aligned  with  the  body  axis  and 
therefore,  must  be  rotated  back  into  the  body  axis  frame.  This  is  accomplished  using 


Feu  —  Pe jj  Le  sin  ($e  -(-  Acre)  i  - (-  Pejj  Le  cos  ( Se  T  Acre)  k 
F eL  F*eL ^e  sin  (6e  T  Ao:e)  i  F*eL Le  cos  ( Se  T  Acre)  k 


(109) 


and 


Fcu  =  PCULC  sin  (Sc  +  Aac)  i  +  PCULC  cos  {5C  +  A ac)  k 
F CL  —  Lc  sin  (<ic  T  Acrc)  i  Lc  cos  ((5C  T  Ao:c)  k 


(110) 


In  terms  of  the  calculation  of  the  moments  produced  by  the  controls,  the  moment  arm  is  also  altered  by  the 
flexible  effects.  For  small  displacements  of  the  beams  relative  to  their  lengths,  the  new  moment  arms  become 


Xcanard  —  ( \J x\anard  (4*/  ( Xcanard  )  T/f,  1 )  ,  &  f  {xCanard)  V f  ,1 


I* elev  — 


^  elev  ^ elev  (-^e  H-  (*^eZeu)  ^7a,l)  ?  %elev  (*^eZeu)  ^7a,l 


(111) 
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VIII.  Results 


The  model  has  been  linearized  at  M  =  8  and  at  an  altitude  of  85,  000  ft  using  velocity  and  flight  path 
angle  as  outputs  and  elevator  deflection  and  total  temperature  addition  in  the  combustor  as  inputs.  One 
point  of  interest  that  can  be  obtained  from  this  simulation  is  the  effect  of  the  unsteady  terms  on  the  poles 
and  zeros  of  the  linearized  system.  Figure  2  shows  a  pole/zero  map  for  the  linearized  rigid  body  system 
when  unsteady  effects  are  not  included,  while  Fig.  3  shows  the  poles  and  zeros  of  the  rigid  body  system  when 
the  unsteady  effects  are  included.  Clearly,  the  unsteady  terms  have  a  significant  effect  on  the  unstable  pole 
and  zero.  Inclusion  of  the  unsteady  terms  makes  the  system  more  unstable,  while  increasing  the  frequency 
of  the  non- minimum  phase  zero.  The  actual  pole  and  zero  locations  are  given  in  Table  1. 


Unsteady  OFF 

Unsteady  OFF 

Unsteady  ON 

Unsteady  ON 

Poles 

Zeros 

Poles 

Zeros 

-0.000137  ±  jO.023358 

1.8694 

—0.000151  ±  j'0. 02332 

2.64 

0.8729 

-1.8694 

1.125 

-2.64 

-0.9159 

0 

-1.216 

0 

-.00107 

N/A 

-.001027 

N/A 

Table  1.  Poles  and  Zeros  of  the  Rigid  Linearized  System. 


Pole-Zero  Map 


Figure  2.  Poles  and  Zeros  of  Linearized  Rigid  Body  System  -  Unsteady  Effects  Off. 


Now,  the  flexible  effects  are  included  in  the  vehicle  model.  Figure  4  shows  a  pole/zero  map  for  the  flexible 
vehicle  when  unsteady  effects  are  not  included.  Figures  5  and  6  show  pole/zero  maps  when  the  unsteady 
effects  are  included.  The  difference  here  is  that  the  stability  derivatives  due  to  ijf.i  and  r]a^  are  included  in 
Fig.  6  while  they  are  not  in  the  results  shown  in  Fig.  5.  In  tabular  form,  the  poles  and  zeros  are  shown  in 
Table  2. 

It  can  be  seen  that  the  unsteady  effects,  computed  using  piston  theory,  significantly  change  the  poles  and 
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Pole-Zero  Map 


Real  Axis 


Figure  3.  Poles  and  Zeros  of  Linearized  Rigid  Body  System  -  Unsteady  Effects  On. 


All  Unsteady  OFF 

Poles 

Zeros 

—  .538  ±  j'26.909 

—  .538  ±  j26.904 

-.36±jl7.99 

—  .36  ±  jl7.99 

1.23 

2.64 

-1.32 

-2.64 

-.00088 

0 

-.00012  ±  j.024 

N/A 

Table  2.  Poles  and  Zeros  of  the  Flexible  Linearized  System  With  All  Unsteady  Effects  Off. 


Unsteady  On  &  Unsteady  Flex.  OFF 

Poles 

Zeros 

—  .538  ±  j26.91 

—.538  ±  j26.91 

—  .36  ±  jl7.99 

—  .36  ±  jl7.99 

1.88 

3.88 

-2.06 

-3.88 

-.00086 

0 

-.00012  ±  j.024 

N/A 

Table  3.  Poles  and  Zeros  of  the  Flexible  Linearized  System  With  Unsteady  Effects  On  and  Flexible  Unsteady 
Effects  Off. 
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Unsteady  On  &  Unsteady  Flex.  ON 

Poles 

Zeros 

—  .55  ±  j'26.91 

—  .56  ±  j26.89 

—  .36  ±  j  17.99 

—  .36  ±  jl7.99 

1.88 

3.88 

-2.06 

-3.88 

-.00086 

0 

-.00012  ±j.024 

N/A 

Table  4.  Poles  and  Zeros  of  the  Flexible  Linearized  System  With  Unsteady  Effects  On  and  Flexible  Unsteady 
Effects  On. 


Pole-Zero  Map 


Real  Axis 


Figure  4.  Poles  and  Zeros  of  Linearized  Flexible  System  -  Unsteady  Effects  Off. 
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Figure  5. 


Figure  6. 


Pole-Zero  Map 


Poles  and  Zeros  of  Linearized  Flexible  System  -  Unsteady  Effects  On,  Unsteady  Flexible  Effects  Off. 


Pole-Zero  Map 


Real  Axis 


Poles  and  Zeros  of  Linearized  Flexible  System  -  Unsteady  Effects  On,  Unsteady  Flexible  Effects  On. 
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zeros  of  the  linearized  system.  Both  the  unstable  pole  and  zero  of  the  linearized  system  move  farther  to  the 
right  in  the  s-plane  when  the  unsteady  effects  are  included.  Only  a  slight  change  in  pole  and  zero  locations 
is  observed  when  the  unsteady  effects  due  to  r\  are  included. 

Figures  7,  8,  and  9  show  the  lift  force,  drag  force,  and  pitching  moment  for  a  one  second  run  of  the 
nonlinear  simulation.  Again,  it  can  be  seen  that  the  unsteady  effects  are  significant  and  are  worth  including 
in  the  nonlinear  simulation. 


x  iQ4  Steady  and  Unsteady  Lift  Force  vs.  Time 


Figure  7.  Lift  Force  (lbs). 


IX.  Conclusions 

In  this  work,  piston  theory  is  used  to  develop  a  model  for  the  longitudinal  dynamics  of  a  2-dimensional 
hypersonic  vehicle  model.  In  particular,  velocities  of  flow  normal  to  the  surface  of  the  vehicle  are  used  in  a 
first  order  piston  theory  framework  to  determine  the  pressures  on  the  surfaces  of  the  vehicle.  The  pressures 
are  then  integrated  over  the  body  to  determine  the  forces  acting  on  the  vehicle.  Piston  theory  is  useful  here 
because  it  allows  the  inclusion  of  unsteady  aerodynamic  effects.  Initial  analysis  has  shown  that  the  unsteady 
effects  are  significant. 
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Steady  and  Unsteady  Drag  Force  vs.  Time 


Figure  8.  Drag  Force  (lbs). 


x  iQ5  Pitching  Moments  vs.  Time 


Figure  9.  Pitching  Moment  (ft-lbs). 
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A.  Flexible  Effects 


From  Eq.  87 
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(112) 


+  AoaC1C3 


„  _  Xa  „  L 
-  sin  Pa,lxa - cos  Pa,lxa - ~ —  ■ 

0a, 1  0?  i 


-  cosh  Pa  ixa  - 

L  ’  t>h 


0a,  1 


0a,  1 

( xa  —  La)  cosh  {^Pa,l  (xa  —  L 


2  AalPcf  . 
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V  J  Pa, 1 


-  La 


sin  (/3a,l  (xa 


l«>)] 

L»>)1 


Aa  l^-'2^-'4  r  _  /  _  \i 

- ^ -  [sin/3ajlxa  -  sin  (0a,i  (®a  -  La))J 


+  A, 


[~xa  1  _  ( xa  -  L 

-  sinh  Pa  i®  a  H - ~ —  cosh/3a  ixa  H - 

0a, 1  Pl^  Pa,  1 


xa  ~  La)  ( 

sinh  I  Pa,l 


Pa,  1 

( xa  ~  La)) 
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where 


Ci  =  sin  Pa,ixa  -  sinh  (3a,\xa 
C2  =  cos  (3a,\x  a  +  cosh/3a>i5a 
P  —  P 

J-  no  J  p. 

C-3  = 


~La 

C4  =  C3  ( La  -  Xa)  -  Pe 


(114) 
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From  Eq.  97 
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From  Eq.  98 
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From  Eq.  99 
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Continued  from  Eq.  99 
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9oo  Sc  V  7  L  ^a.l  f^,l  J 

—  1  /  q  3  _  2  _2  __  _ 

— —  \Pa  1  xa  sin^1  0a,lxa  ~  3Pa  1  xa  cosh  0a,lxa  H-  ®0a,lxa  0a,lxa  ~  6  cosh  0a,lXa 

a,  1 


P.aacosjrjn+Tj) 

- - - ~Aa,l  (cos  0a,lxa  ~  cosh /3 a,lXa) 


Qoo  Sc 

[pi  1  (*a  -  L„)3  sinh  /3 

^,1  ' 

cos  (ti  [/  +  T2) 

- ; - ^ - J5Aa,l 


goo  Sc 


'o ,  1  (*a  -  La)  -  302a  l  (. 

0a,  lx, 


paaa  cos  (tx  tj  +  r2) 

- - - ~Aa  1  (  c°s  /3a  jxa  —  cosh /3a  ^5 

goo  Sc 

xa  —  La)2  cosh  0a>i  ( xa  —  La )  +  6/3a  ^  (Sa  —  La )  sinh  (Sa  —  J^a)  —  6  cosh  Sa,l  (®a  —  ^ 


25a  _  0a  lxl  ~  2  _  f  2  (*a  -  La)  _  0a  l  (xa  -  La)2  -  2  V 

— ^ —  sin  0a,1xa - ’—s -  cos  0a  \xa  ~  {  ^ -  sin  0a  \  ( xa  -  La) - 1 - ^ -  cos  0a  i  ( xa  -  La)  > 

^a.l  Pi.  1  l  ^.1  ^a.l  J. 


cos  (rlj!7  +  ra) 

+ - — - J5Aa,l* 

Qoo  Sc 

<  [^S2  i  *2  cosh  0a,  i  Sa  -  2/3a  ^  5  a  sinh  /3a,  1  xa  +  2  cosh  1  xa  ~  |i^a  1  (*a  —  cosh  ^  (5a  —  i^a)  —  2Sa,  1  (xa  ~  La)  sinh  0a  ^  (5a  —  La )  +  2  cosh  /3a>  4  (Sa  — 

(3?  2  f  2(xa  -  La)cos/3a  !  (x0  -  Z,a)  ^  i  (*«  “  La  )2  -  2 

- ^ - ! - 1 - ‘ - 5 -  sin  0a ,  1  ( xa  -  La) 

32  ,  0$  ± 


cos  (ti  n-  +  t2)  25a  cos  /3_  iia  ^a  l^a  _  2 

1 - 1 - J5Aa,l  - ^ - 1 - 1 - -  sin /?a,  i®  a  —  {  ~ 

g°o  Sc  ^a,l  03, 1  l 


cos  ('rl,C7  +  T2)  r  2  _2  _  _  _  _  1 

- J5  Aa  I  0a  1xa  sinh  0a  \xa  ~  '20a,lxa  cosh  0a,lxa  +  2  sinh  0a  \xa\ 

QooSc  ’ 

cos  (Ti  n  +  r2 )  r  2  9  i 

H - - - —  J$Aa,  1  f^a  1  _  La)  sinh  0a,l  (xa  ~  La)  ~  20a,l  (xa  ~  La)  cosh  0a,l  (xa  ~  La)  +  2  sinh  0a,l  (xa  ~  La)\ 

gooSc  ’ 

cos  (Ti,u  +  T2)  r  _  _  r  _  _  n 

\ - - - - - -  ^6^a,  1  I2  -  cos  0a,lLa  +  0 a, 1  (xa  ~  La)  sin  0a,lLa  +  0 a, 1  ixa  -  La)  |sin  0a^\xa  cosh  0a,l  (xa  ~  La)  +  cos  0a^\xa  sinh  0a,l  (xa  ~  La)> J 

QooSc 

\xa  sinh  0a,  1  (xa  ~~  La)  —  cos  0a,lxa  cosh  0a,l  (^a  —  La)  +  sin  0a,l  (^a  —  La)  sinh  0a,l^aj 


cos  (t1jC7  +  T2) 

H - - - —  J§Aa  4  2  cos  0a  jia  cosh  0a  \xa  —  sin  0a  \xa 

Qoo  Sc 

cos  (Ti  r;  +  i-2)  r 

H - - - —  JcAr,  i  —  cos  0a  i  (xa  —  La )  cosh  /3a  iia  —  0a  1 

QooSc  ’  L  ’  ’ 

COS  (ti  tj  +  r2) 

H - - - ;-JgAa  —  0a  \  ( xa  —  La)  sin  0a, 1  (xa  ~  La  )  cosh  0a  ^  xa  —  0a  ^  (®a  —  La  )  sinh  0a>  i  La  —  cosh  0a  ^ L 

QooSc  L 

(119) 


(xa  —  La)  cos  0a  ^  (xa  —  La)  sinh  0a,lxaj 


where 


acw  =  |  dcM 

^Va,  1  1  ^Va,  1  p2 


(120) 


•/i  =  paaami 

J2  =  -2kpaaami  +  paaa  {mi5a  +  m2}  +  (peaa  +  paae)  mi 
d3  =  paaamik 2  -  2 kpaaa  {miX0  +  m2}  -  {peaa  +  paae}  rn xk  +  {peaa  +  paae}  {m\Xa  +  m2}  +  /9eaem.i 
Ji  =  Pa&a  {mixa  +  m2}  k2  -  {peaa  +  paae}  {miia  +  m2}  k  +  peae  {mixa  +  ?n2} 

t/5  2kpaaa  [Pe^a  ~k  Pa^e) 

J&  =  PaCLak 2  -  k  (peCLaPaCle)  +  peO,e 

_  Poo-Pe 

ka  —La 

n  —  aoo  ae 

a  ~  ~La 
h  —  La  xa 
mi  =  tan  (n^  +  r2) 
to2  =  —  itanTi  u 


(121) 
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